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Abstract 

Ratios of integrals can be bounded in terms of ratios of integrands under certain mono¬ 
tonicity conditions. This result, related with L’Hopital’s monotone rule, can be used to 
obtain sharp bounds for cumulative distribution functions. We consider the case of non¬ 
central cumulative gamma and beta distributions. Three different types of sharp bounds 
for the noncentral gamma distributions (also called Marcum functions) are obtained in 
terms of modified Bessel functions and one additional type of function: a second modified 
Bessel function, two error functions or one incomplete gamma function. For the noncen¬ 
tral beta case the bounds are expressed in terms of Kummer functions and one additional 
Kummer function or an incomplete beta function. These bounds improve previous results 
with respect to their range of application and/or its sharpness. 

2000MSC: 33E20,33B20,26D07,26D15 

1 Introduction 

Given two functions fi{x) integrable in [a, 6] and continuous in (a, &) and 

= f fi{t)dt, i = 1,2, (1) 

J a 

we have that, on account of L’Hopital’s rule, the ratio F 2 {x)/Fi{x) has the same limit as 
f 2 {x)/fi{x) for cc —>■ 0 + and the same is true for (E 2 (x) — F{b))/{Fi{x) — F{b)) as x —>■ b~. 
Additionally, as we will see, when / 2 (x)//i(x) is continuous and monotonic, the ratio becomes 
a bound for the ratio of integrals. This result is closely related to L’Hopital monotone rule [5J 
Theorem 2] (see also [1] Lemma 2.2] and [6l p. 42]). 

These bounds on ratios of integrals can be used for bounding cumulative distribution func¬ 
tions Gi{x) = Fi{x)/Fi{b) when additional information is available, like, for instance, the values 
of the difference G 2 (x) — Gi(x) or a known explicit expression for one of the functions Gi{x). 
We consider this type of bounds for the noncentral gamma and beta distributions, obtaining 
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different families of sharp bounds for different selections of the pair of functions Gi. The non¬ 
central gamma distributions are also called Marcum functions and have as a particular case 
(central case) the incomplete gamma function ratios. On the other hand, the central case for 
the beta distribution is the incomplete beta function ratio. 

The gamma and beta distributions are classical cumulative distributions appearing in many 
scientific applications. The noncentral gamma distribution is an important function in radar 
detection and communications and it is widely used in statistics and probability theory (see for 
instance m and references cited therein). Beta distributions are also important functions in 
statistics and probability. The computation of these functions, particularly for the noncentral 
case, is difficult. Only recently accurate algorithms for both the lower and upper tail noncentral 
gamma distributions (Marcum functions P and Q), valid for a large range of parameters, were 
developed. For the noncentral beta case the situation is worse and the available methods of 
computation are based on the definition of the lower tail distribution in terms of incomplete 
gamma functions (Eq. (jl6|) b Probably due to this difficulty, many researchers have been 
involved in obtaining bounds and approximations in terms of simpler functions, particularly in 
the case of the gamma distribution (see, for example, [3] and references cited therein). 

We obtain several types of bounds for the gamma and beta distributions. For the noncentral 
gamma distribution we obtain three types of bounds in terms of modified Bessel functions alone, 
in terms of a modified Bessel function plus two error functions or in terms of a Bessel function 
and an incomplete gamma function. The first type of bounds was recently considered in m 
but using a different approach, while the second type is related to some of the bounds in [3]. 
For the noncentral beta case the bounds are expressed in terms of Kummer functions alone or 
in terms of an incomplete beta function and a Kummer function. In addition to the bounds 
related to L’Hopital rule, we will see how the recurrence relations satisfied by these functions 
provide complementary bounds. Both for the gamma and beta cases, the bounds we obtain 
improve previous results with respect to their range of application and/or its sharpness. 

We will also prove some monotonicity properties and bounds for Kummer functions which 
are needed in the construction of the bounds for the beta function (see the appendix). Ad¬ 
ditionally, we will briefly explain how the bounds are useful for computing the inverse of the 
distributions (named quantile functions in statistics). 


2 Bounds related to L’Hopital’s rule 

In this section we consider functions defined by integrals. The same results are valid in general 
for functions Fi{x) satisfying F)(a) = 0 and which are continuous in [a,b] and differentiable in 
(a, 6); we are not loosing generality by writing these functions as integrals. We assume that 
Fi{b) ^ 0 and we will deal with normalized integrals Gi{x) = Fi{x)/Fi{b) (so that Gi{b) = 
G 2 (b) = I); it is clear how to write the results for Fi from those for Gi. 


2.1 Bounds for ratios of integrals 

The following theorem can proved by using Rohe’s theorem and Cauchy’s mean value theorem. 


Theorem 1 Let 


Gi{x) = / gi{t)dt, I = 1,2 
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with Gi{b) = G 2 {b) = 1, gi(x) and g 2 {x) integrable in [a, 6] and continuous in (a,b) and with 
gi{x) ^ 0 in (a,b). Let Gi{x) = 1 — Gi(x). 

Then, if g 2 {x)/gi(x) is strictly monotonic in (a,b) the following holds: 

1. There exists one and only one xq in (a,b) such that gi(x(f) = g 2 {x(f) 

2. Leth(x) = G 2 {x)/G i{x) andr{x) = g 2 {x)/gi(x). The following statements are equivalent 
in {a, b): 

(a) h{x) < 1 (h{x) > \) 

(b) r(x) is inereasing (decreasing). 

(c) h{x) < r{x) (h{x) > r{x)) 

(d) h{x) is increasing (decreasing) 

The same holds for h{x) = G 2 (x)/Gi(x) but reversing the inequalities in (a) and (c). 

Proof. First we notice that Gi{x) are continuous functions in [a, &] and differentiable in 
{a,b). 

1. Let D{x) = G 2 {x) — Gi{x). Because D{a) = D{b) = 0, by Rolle’s theorem there exists 
xo G {a,b) such that D'{xo) = 0, and therefore g 2 ixo) = gi(xo). And because g 2 {x)/gi(x) is 
strictly monotonic no other value x ^ xq exists such that g 2 (x)/gi(x) = 1. 

2. We prove the chain of implications (a) ^ (6) (c) (d) (a) for the case h{x) < 1; 

the case h{x) > 1 is analogous. 


1. (a) (b): We apply Cauchy’s mean value theorem for Gi(x), i = 1,2 in [a,a;o], with xq 

the only value such that r(xo) = 1- Then there exists c G (a,xo) such that 

Gi{x) G,{x)-G,ia) gi(c) 

and because h{x) < 1 then r(c) < 1 with c < xq- And because r{xo) = 1 and r(x) is 
strictly monotonic by hypothesis, then it must be strictly increasing. 


2 . 


(b) => (c): We consider Cauchy’s mean value theorem in [a, x], x < b: there exist c G (a, x) 
such that 


h{x) 


G2{x) 

Gi{x) 


92{c) 

9i{c) 


92{x) 

9i{x) 


= r{x), 


( 2 ) 


where the last inequality holds because c < x and r{x) is strictly increasing. 


3. (c) (d): Taking the derivative 


h'{x) = {r{x) - h{x)). (3) 

Gi[x) 

And because h{x) < r{x) then h{x) has positive derivative (because gi{x) ^ 0 for any 
X G (a, 6) then gi{x) and Gi(a;) have the same sign). 

4. (d) (a): Because h{x) is strictly increasing in (a, 5) and h{b) = 1 then necessarily 

h{x) < 1 in (a, b). 
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For Gi{x) we can proceed analogously. For instance, for the second implication we can 
apply Cauchy’s mean value theorem for Gi{x) in [x, h\. 

G2{b) - G2 {x) _ G2{x) _ 32(c) 32 ( 3 ^) 

Gi{b)-G 2 {x) Gi{x) 31 (c) 31 (a;)’ 

The third implication follows by considering the derivative of h{x) = G' 2 (x)/G'i(x), which is 

Hx) - h{x)) . 

□ 


Remark 1 We have proved that g 2 {x)/gi{x) is a bound. In addition, there is the trivial bound 
1 (because r(x) < 1 or r{x) > I). Therefore if, for instance, h(x) <\ we have 


G2{x) 

Glix) 


< min 


, 92{x) \ 

’ 51 ( 3 :)/ 


G2ix) 

Glix) 


> max 


, 92ix) \ 
’51(3:)/ 


In all cases, the bound 32 (x) /31 (x) is sharp for 6*2 (x) /Gi (x) as x —>■ 0 + and /or G 2 (x) /Gi (x) 
as X ^ b~. The trivial bound 1 is sharp on the opposite extreme of the interval. The bound 
32 (x)/ 3 i(x) for G2(x)/Gi(x) ((52 (x)/Gi (x)^ is better than the trivial bound only for x < xq 
( x > Xo). 


Remark 2 That the monotonicity properties of 32 (x)/ 3 i(x) are inherited by G 2 (x)/Gi(x) 
(implication ( 6 ) (d )) is the result named monotone L ’Hopital rule Theorem 2], and it is 
proved in Lemma 2.2] similarly as we have proved the two implications (b) => (c) (d). 

Remark 3 An alternative proof for (a) => (b) => (c) => (d) can be given which provides a 
graphical picture of the situation. The proof goes as follows: 

We start from (a) and prove (b) as before. Now, we use ^ to conclude that h'(x) has 
the same sign as r(x) — h(x). With this it is easy to check that the graph of h{x) is inside 
the region {a,b) x (h(a'^),l) and below the graph ofr{x). Indeed h(b~) < 1 < r{b~) (because 
h{x) = G 2 (x)/Gi(x) < 1 and r{b~) > r(xo) = Ij and then, considering ( 0 ), h(x) is increasing 
close to b; but it stays being increasing in the rest of the interval because, as we move to the 
left of the interval, the graph of h{x) can not cross the graph of the increasing function r(x). 
If it did at some x = Xc < xq, then h(x) would become decreasing for x < Xc, implying that 
h{x) > h(xc) = r{xc). Because Xc > a and r(x) is increasing, this would contradict the fact 
that h{a'^) = r(a'^) on account of L’Hopital’s rule. 


2.2 Bounds for the functions 

The bounds for the ratios can be translated into bounds for the functions if some of the 
functions Gi or G 2 are explicitly know. We will use this case for obtaining bounds relating 
the noncentral and the central distributions; also, for the gamma distribution we will obtain 
bounds in terms of error functions using this idea. 

Another possibility is that there is some relation between Gi and G 2 as, for example, the 
difference Gi(x) — G 2 (x). This happens both for the gamma and beta distributions, where the 
recurrence relation they satisfy provides such a relation. 
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Considering then that Gi{x) — G 2 {x) is known we write: 


Gi{x) 


G^{x)-G2{x) 
. G2{x) 

Gi{x) 


and considering that the case h{x) < 1 of theorem [2] holds we have G 2 {x)/Gi{x) < r{x) = 
92 {x)/gi{x), with g 2 {x)/gi{x) < 1 if x < ccq ( 92 ( 2 : 0 )/ffi(a:o) = !)• Therefore 


Gi(x) < 


Gi{x)-G2{x) 
1 — r{x) 


= C/i(x), X < xo, 


(4) 


and similarly 

G 2 < U 2 {x) = r(x)C/i(x), X < xq. (5) 

These two bounds are sharp as x —>■ a+. 

In a similar way, we obtain 


Gi(x) < —17i(x), X > xo, 


( 6 ) 


and 


G2(x) < —U2{x), X > Xo, 


(7) 


and these two bounds are sharp as x —>■ 6 . 

Next we consider the application of Theorem [T] to the gamma and beta distributions. 


3 Bounds for gamma distributions 

The Marcum functions (or, equivalently, non-central gamma or chi-squared distributions) are 
the cumulative distribution functions 

ry 

Py.{x,y)= / 9 ^(x,t)dt, Q^(x,y) = 1 - P^(x,y) 

20 

with 

9ii{x,y)=[£^ e %_i(2y^). 

We have P^{x, 0) = 1 and P^{x, -l-oo) = 1, as corresponds to a cumulative distribution function 
in [0, -boo). 

Remark 4 We observe that the role of the variable x in section\^is now played by the variable 
y. We use this notation in order to be consistent with previous notation. This notation is used 
in the rest of this article. 

Integrating by parts we obtain the recursion formulas: 

Pu+ii^^y) = Pvi^^y) - 9u+ii^^y)^ Qu+^i^^y) = Qni^^y) + 9u+ii^^y) ( 8 ) 

We consider three types of bounds. The first type of bounds were already explored [15] (but 
using a different approach); some of the results in that paper are particular cases of these type 
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of bounds. These bounds involve monotonicity properties over /i and a know expression for the 
difference given by the recurrence relation. The second type of bounds (error function bounds) 
can be obtained from the monotonicity over ^ and the particular expression for /r = 1 / 2 , 
/r = 3/2. Finally the third type of bounds (in terms of incomplete gamma functions) are based 
on the monotonicity with respect to x and the particular value for a; = 0 . 


3.1 Bounds for the functions nsing their recnrrence 

The hypotheses of Theorem 12.11 (case h{x) < 1) hold with the correspondences Gi —^ 

(Gi Qfj), gi gfj,, G 2 -)> P^^+i (G 2 -)> Qf_t+i) and 32 gti+i An because g^+i(a:, y)/gf^{x, y) 
is monotonic as a function of y [m lemma 2 ], we have: 


P 




'x,y) 


Pfj^ix,y) 


< 


and 


gf,+i{x,y) ^ 
9ti.i.x,y) 

Qii.+i{x,y) 

Qt^{x,y) 




X If_y-i{2^/xy) 


> Cf,[x,y). 


The first bound is slightly improved in Theorem 3 of m while the second is given in Theorem 
6 of that reference. 

From these bounds for the ratios, we can obtain bounds for the function considering the 
expression given for Gi — G 2 which is obtained from the recurrences ([HI) . Therefore the bounds 
(IHTfl) hold, with the correspondence Gi(a:) —G 2 (a:) —>■ 5^+1 (x, y). These are some of the bounds 
in section 3.2 of m- 

From the construction of the bounds (related to L’Hopital’s rule) we deduce they are sharp 
as 2 / — >■ 0 for P^ and as ?/ —>■ +00 for In addition, as explained in m, the bound for P^ is 
also sharp as x —)> +00 and y + 00 . 

Additional bounds can be obtained from the recurrences alone (| 8 |), without invoking The¬ 
orem [T] Indeed, because the P and the Q functions are positive P^(x, y) > y^+i(x, y) and 
Qfj.(,x,y) > g^(x,y). We can iterate the recurrence and obtain more sharp bounds. For in¬ 
stance P^{x,y) > y) We refer to [15] for further details. 


3.2 Bounds in terms of error functions 


We take G 2 (y) = g^+a{x,t)dt and Gi(y) = g^{x,t)dt, then 


r(y) = 


9fj.+aix,y) _ , , ,Q, /^_i+Q,(2yGiy) 


9i^ix,y) 


= C{a,x){^y 


t/i —1 (2\/ xy) 


with G(a, x) not depending on y. The ratio r(y) is strictly monotonic (increasing) as a function 
of y for y > 1 as the next lemma shows. 

Lemma 1 The function is increasing as a function of t for a, y > 0. 

Proof Using the relations in (TUI 10.29(i)] we obtain: 






(t)\ f lu+a-lit) lu-iii') 






it) 


I nit) 
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And because /I^{t) is increasing as a function of > 0 [TH Lemma 2] the result is 

proven. □ 

Now, using that r(j/) is increasing and applying theorem 12.II we have that 


Theorem 2 If v > ^jl > 1 then 


and 


Qu{.x,y) > 


I^-i{2y/xy) 

/^_i(2yiy) 

I^-i{2^/3yy) 

Ifj,_i{2y/¥y) 


Pt,{x,y) 

Qp.{.x,y) 


Remark 5 We conjecture that t > 0 is increasing as a function of y > —1/2 

and, therefore, that lemma\l\ holds for y > —1/2, a > 0 and theorem\^ holds for fj, > 1/2. A 
numerical check shows that this is true. 


Taking /r = 1/2 in the previous theorem (and taking into account the previous remark) we 
obtain bounds similar to those in [J. For instance, and in terms of the variables a = 
b = (as in [3]). From the bound for Q^, in the previous theorem we have: 


Qi,(a, b) > b) 


B]}\a,b) = 


and from the bound for Pi, 


G,(c 


cosh(a6) 


i (w) + (^)) 


(9) 


Qy{a,b) > B^v\a,b) 

_ 1 _ Gt/(a, b) 


B\f\a,b) = 1- 


cosh(a6) 


(erf( 


b—a \ 

V 2 ) 


erf 


(^)) 


i{ab) 


( 10 ) 


where Q„{a,b) = Q{af 12,1? 12') (Q,^(a,b) is the function denoted as Q„(a,b) in m) 


GAa,b) 


a jfv 
2\2\a) 


The bounds are valid for ^ > 1/2 and any a > 0, 6 > 0. 

As a marginal note, we notice that the inequality ([9l) implies the following relation between 
error functions and Bessel functions: 


('* (V)+(^)) ^ 

valid for v > —1/2, a,b > 0. 

The bounds and m can be compared with the results is [3] (eqs. (8) and (17)). For 
instance, Eq. (8) of the reference is: 


Q 1 /(u, b) ^ 


G.(q,&) 

sinh(a6) 


erfc 


b — a 


— erfc 


72 


4_i(a&), 
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valid for 6 > a > 0 and v > 1. Which is slightly superior to © when it holds, but of more 
restricted validity; only when a is small a < 2 and b is small (but b > a) the advantage becomes 
noticeable. On the other hand, in the comparison between (fTUl) and Eq. (17) of [3], we observe 
that is a better bound except if a < 6 is quite close to b. 

We can obtain sharper bounds (but of more restricted validity with respect to /i) by taking 
larger fi in Theorem [3] For instance, taking /r = 3/2 we obtain bounds valid for ^ > 3/2 and 
any a,b > 0 which are sharper that the bounds of eqs. (8) and (17) of [3]. Using Theorem [2] 
together with the recurrence we have 


and 


Qb'^ ^ 
Quia, b) > I 


aBlHa,b) ^ 
5tanh(a6) 

(a, b) 
6 tanh(a6) 




valid for i/ > 3/2. 


( 11 ) 

( 12 ) 


3.3 Bounds in terms of incomplete gamma fnnctions 

We take G 2 (j/) = gvip^x,t)dt and Gi{y) = {x, t)dt, with p > 0, p ^ 1, then 


riy) 


9uip^x,y) 

9uix,y) 


G.(p, 


Iu-l{2py/^) 

Iu-ii2y/xy) 


where Cuip,x) does not depend on y. Therefore, riy) has the same monotonicity properties 
as /^(pz)//^(z) as a function of z. Using [TOl 10.29(i)] we obtain 

f It^ipz) \ ^ Ifiipz) f 1^,+iipz) _ If,+i{z) \ 

V V If^ipz) If^iz) ) ’ 

and because xl^+iix)/1^{x) is increasing as a function of x for p > — 1 (see for instance [M] 
Lemma 2]), riy) is increasing for > 0 if p > 1 and decreasing if p < 1. Applying theorem 12.II 
for p < 1 we obtain 



Puix,y) < y-^y^.Puip^x,y) 

Iu-ii2py/xy) 

(13) 

and 

Quix,y) > pM-lg-.d-P^) Q,(p2a;, y) 

Iu-ii2py/xy) 

(14) 

and taking the limit p —>■ 0 



Puix,y)<e ^ (y/xy)^ Iu-ii2y/xy)juiy) = b\}\x,y), 
Quix,y) > e““ (\/5p)^ " Iu-ii2y/¥y)ruiy) = b^u\x,y) 

(15) 


where 












are the lower and upper incomplete gamma functions. 

Observe that, as expected, 

lim 6W(a:,2/) = = P^{0,y), 

a:— ^ ) 

lim 6[?l(a:,2/) = :p^r^(?/) =Q^{0,y) 

X-S-0+ 1 [I') 

Furthermore, because Q^(x,y) < 1, the second inequality of (fT51) implies the curious in¬ 
equality 

h{2y/¥y)T^+i{y) < {y/xyY ,v > -l,a:, y > 0 

Bounds in terms of incomplete gamma functions are also given in |12j . Those bounds can 
be applied for a = or b — y/2y not much larger than 1. Our bounds do not have this 
limitation: The bounds |12] appear to sharper for small a or & but our bounds are sharper for 
larger parameters. 

4 Bounds for beta distributions 

The cumulative noncentral beta distribution can be defined as 

OO - 

Ba,b{x, y) = ^ ^ (ly Iy{a + j, b) (16) 

where /y(a, h) is the central beta distribution 

= = (17) 

B[a,b) Jo T[a + b) 

A different notation used in the literature (see for instance 0) is Ix{a, b, A), where Ba,b{x, y) = 
Iy{a, b, x). The complementary function is defined by Ba,b{x, y) = 1 — Ba,b{x, y). 

Using (HU in (Hii) we get 

ry 

Ba^h{x,y)= / ga,h{x,t)dt (18) 

with 

_ 31 ^ 2 

5a.b(a:, y) = (a + b, a, y) (19) 

and similarly for Ba,b{x,y) but with the integration going from y to 1. 

As corresponds to a cumulative distribution function, we have 

0 = Ba,bix, 0 ) < Ba,b{x, y) < BaY^, 1 ) = 1 

for a, & > 0 and x,y > 0. Notice that, as for the gamma distribution, the variable y corresponds 
to the variable x of section [2l also, as for the gamma distribution, x is called the noncentrality 
parameter. The central case corresponds to a: = 0. 
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4.1 Bounds for the functions using their recurrence 

Using the recurrence relations for the central distribution given in [TT] 8.17(iv)] together with 
the definition (ITCl) it is easy to get analogous relations for the noncentral case. 

Using [TT] 8.17.20] we obtain 

Ba,b{x,y) = Ba+i,b{x,y) + + 5, a + l,xyl2), (20) 

a 

from [IS 8.17.21] 

Ba,bix,y) = Ba,b+iix,y) - M(a + b, a,xyl2), (21) 

0 

and from [TT] 8.17.19] 

Ba,bix,y) = Ba-i,b+i{x,y) - jVf(q _|_ b,a,xy/2), (22) 

by 

where 

For the complementary function, Bafi{x,y) = 1 — Bafi{x,y) the same recurrences hold but 
with the sign of the inhomogeneous term reversed. From these recurrences, it becomes clear 
how the results of section [5] can be applied in this case. We give the inequalities for ratios of B- 
functions; for the S-function the same holds but with the inequality reversed. We concentrate 
on the recurrence relation (|20l) : for the rest of recurrences, a similar analysis is possible. 

From (EiH) we observe that the role of Gi and G 2 in section[2]are played by Ba^b and i?a+i.b- 
Therefore 

Ba+i,bix,y) ga+i^b{x,y) ^ g + b M{a + b+l,a+l,z) 

Ba,b{x,y) ga,b{x,y) “ ^ M{a + b,a,z) '' 

where z = xyj2, provided that the function ratio on the right-hand side is monotonic, which is 
obviously true for the central case a; = 0. The proof of monotonicity can be obtained by using 
similar ideas as those used in m for the ratios of modified Bessel functions. In the Appendix 
we prove that zM{a + b+ 1, a + 1, z)/M{a + b,a, z) is increasing; in addition, using the bound 
given in this theorem, we obtain the simpler bound 

B^i^b{x,y) .^(2; + 1 _ a)2 _|_ 4(a _|_ b)z^ , 2; = xy /2 

BaAx.y) X \ J 

Now we obtain bounds of incomplete beta functions from the bounds for the ratios as 
described in section [521 Using (1501) and ((55)) we have the bound for Gi 


and for G2 


Ba,b{x,y) < 


Ca,b{x,y)M{a + b,a+ l,xy/2) 


1 - 


9a+l,b 
9a,b 


Ba+iMix,y) < 


Gg^bjx, y)M{a + b,a+l, xy/2) 

( 9a,b 

a -^- 

\5a-|-l,h 


(25) 


(26) 
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Let us particularize these bounds for the central case (a; = 0). The inequality (E5|) gives 


Ba,bi^,y) < 




y < 


B{a,b){a — {a + b)y)’ a + b 


and (1^ 


BaA^^y) < 




,y < 


a — 1 


B(a, b){a — 1 — (a — 1 + b)y) ’ a — 1 + 5 

which is essentially the same bound. 

For Ba,b{0,y), we have similarly to (E7)) 


Ba,bi0,y) < 


vA^-yf 


y > 


B{a,b){{ 0 ‘+ b)y — a)^ a + b 


( 27 ) 


(28) 


(29) 


and the bound related with (1^ . 

The bound (1251) compares favourably with the upper bound in [4l Thm. 5]. It is simpler 
and sharper except in some cases when a is small (a < 3). For large values of a and/or b our 
bound is several orders of magnitude sharper. 


4.1.1 Bounds using only the recurrence 

The recurrence relations, without invoking Theorem [TJ provide sharp bounds themselves. For 
instance, considering (l20l) we conclude that 


Ba,bA^y) = ^W (a + 6, o + 1, xy/2) 

a 

and for the complementary function we have 

Ba,bA^y) > M(a + b- I,a,xy/2). 

a — 1 

For the case of Ba.b(x^y) we can construct a convergent sequence of lower bounds by 
applying the recurrences consecutively. In that way, after applying N times the recurrence we 
have 


Ba.bA^y) = Ba+N+i,b + y] y^ M (a + b + j,a+ j + l,xy/2) 


3=0 


(a + l)j 


which gives the bound 


Ba,b{x,y) > 


CaAx,y) 

a 


E 


3=0 


(a + ^)j 

(a + l)j 


y^M(a + b + j,a + j + l,xy/2) 


(30) 
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4.1.2 An application: quantile function estimation for the central beta case 

The previous bounds can be used as estimation for the cumulative distribution, particularly 
when they are sharper (at the tails). But they can also be used for computing estimations of 
the inverse distribution or quantile function. For instance, the bounds for the central gamma 
distribution were used in some cases as starting values for the inversion method described in 
[8]. Similarly, the bounds for the central beta distribution are useful as starting values for 
inverting the equation Ba^b{^,y) = P, with f] small. 

Considering the estimation (EZl) and because Bafi{0,y) is increasing as a function of y and 


f{y) B{a,b){a — {a + b)y) 

is also increasing as a function oiy iiy < a/(a+b), then the solution of f{y) = /3 in (0, a/(a+b)) 
gives a value yi such that Ba,b{0,yi) < P- Therefore the value yp such that Ba,b{0,yp) = /3 is 
such that yi < yp and it should be close to yi if /3 is sufficient small. Depending on the values 
of a and 6, the lower bound yi guarantees convergence of the fourth order Schwartzian-Newton 
method [13] to yp; and when it does not guarantee convergence, an upper bound will do. 

An upper bound for yp can be obtained from the lower bound for Ba^ix^y) (1501) . For 
instance, considering N = \ (and a; = 0) we have 


Ba,bP^,y) > 


aB{a, b) 



and the solution of g{y) = p gives an upper bound for yp, yu > yp- Therefore yp S (?/;,?/„). 

The equations f{y) = p or g{y) = a can be solved by any numerical method for nonlinear 
equations. A possibility with good convergence properties is to use the fixed point method 
obtained by isolating y as follows 

y = h{y) = iPBia, b){a - (a + b)y){l - 


we can iterate j/n+i = h{yn) starting from j/o = 0 in order to solve numerically the equation 
f{y) = p. We can proceed similarly for g{y) = P taking now 

Ky) = (^ + {aPB{a, 6)(1 - 

These approximations are used as starting values for the the inversion of the cumulative 
beta distribution for small values of ,5 |7]. One of them (either yi or j/„) guarantees convergence 
of the fourth order Schwartzian-Newton method to yp. 


4.2 Bounds in terms of the central distribntion 

The roles of Gi and G 2 are now played by ga,b{px, y) and ga,b{x, y). We have 


ga,b{px,y) 

9a,b{x,y) 


C{a, b, x) 


M{a + b, a, pz) 
M{a + b, a, z) 


z = xy/2 
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and C{a, b, x) does not depend on y. Therefore the monotonicity properties of r{y) = ga,b{px, y)/ga,b{x, y) 
as a function of y are the same as the properties as a function of z of q{z) = M {a, /3, pz) /M {a, /3, z), 
a = a + b, j5 = a. Now, we the aid of [21 13.3] we obtain 


aM(a,/3,pz) M(a + l,l3,pz) M{a + l,f},z) 
z M{a,l3,z) _ M{a, P, pz) M{a,/3,z) 


In the appendix we will prove that increasing as a function of z for 

a, /3 > 0. Therefore r[y) is decreasing as a function of y for p < 1 and we have 


Ba,b{px,y) ^ ga^b{px,y) ^ ^_^(p_-^y2 M{a + b,a,pxy/2) 
Ba,b{x,y) ga,b{x,y) M{a + b,a,xy/2) 

An particularizing for p = 0, we have 


Ba,b{x,y) <e ^/^M{a + b,a,xy/2)Ba,biO,y) 


Similarly 

Ba,b{x,y) > e~'^^^M{a + b,a,xy/2)Ba,biO,y), 

which implies that 

5a.&(0,i/)M(a + 6,a,xp/2) < 

for a,b,x > 0, y G (0,1). 


Appendix 

In order to obtain a bound for the beta distributions, we used some monotonicity properties 
for the ratios of Kummer functions that we next prove. 

First, we used the property that zM{a +!,&+!, z)/M{a, b, z) is increasing as a function 
of z. Next we prove this result and we supplement it with additional information. 


Theorem 3 The function 


h{z) 


Af ((2-t-l,&Tl,z) 
M{a,b,z) ’ 


is decreasing as a function of z > 0 if a > b > 0, constant if a 
0 < a < b, while zh(z) is increasing for all a,b > 0. 

If a > b > 0 the following holds 


b > 0 and increasing if 


D{a, b, z) < yh{z) < D{a, b — l,z) 
b 


where 

D{a, b,z) = ^ (^z - b + \/{z - b)"^ + 4azj . 

The upper bound also holds for 0 < a < b, but a sharper bound is given by 


^h{z) < D{a,b,z) 
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Proof. We start with the function zh{z) = ^H{z) with 

M'{a,b,z) 


H{z) = z- 


M{a, b, z) 


taking the derivative and using the Kummer equation zM''{a,b, z) + {b — z)M'{a,b, z) — 
aM{a, b,z) = 0 we find 


H' (yz) — a + [ 1 + 


1-6 


H{z) - -H{zf 
z 


We compute the characteristic roots A(z) (such that if H{z) = X{z) then H'{z) = 0) by solving 
a + ^ ^ A(z) — ^A(z)^ = 0. Because o > 0 one solution is positive and the other 

negative. The interesting solution is the positive, which is A+(z) = zD{a, 6—1, z). 

Because a, 6 > 0, H{z) > 0. We will have H'{z) > 0 when 0 < H{z) < A+(z) and H'{z) < 0 
when H{z) > A+( 2 ;), but the second possibility can be easily ruled out. 

We have A+(0+) = 0+ and A(^(z) > 0, z > 0. On the other hand, i?(0+) = 0+ and 
i/'(0+) = 1, which implies that i?(0+) < A+(0), but then necessarily H{z) < A+(z) and 
H'{z) > 0 for all z > 0. Indeed, the graph of H{z) can not intersect the graph of the 
characteristic root A+(z) because A+(z) is increasing and i?(0+) < A+(0). 

This proves the monotonicity properties of zh{z) and the upper bound of the theorem. 


Now we consider G(z) = 


M'{a, 6 , z) 


equation 


M (a, 6 , z) 


= ^h{z). Taking the derivative and using Rummer’s 


G'(z) = -+ 1-- G(z)-G(z) 


For a > 0 there is only one positive characteristic root, which is A+(z) = D{a,b,z), which 
is decreasing if a > 6 > 0 , equal to 1 if a = 6 and increasing if 0 < a < 6 ; in addition 
A+(0+) = a/ 6 . On the other hand G(0+) = a /6 with G'(0+) < 0 if a > 6 , G'(z) = 0 if a = 6 
and G'(0'^) > 0 if a < 6 . With this, and using similar qualitative arguments as before, the 
monotonicity properties for h{z) and bound is proved (which gives the lower bound in the 
theorem). □ 

For obtaining the bounds for small x, we used monotonicity of the ratio M (a+1, b,x)/M (a, 6 , x). 
We prove this result and an associated bound: 


Theorem 4 The function 

M(a + 1 , 6 , z) 

M(a, 6 , z) 

is increasing as a function of z for a, 6 , z > 0 and the following bound holds for that range of 
the variables: 


M{a + 1, 6 , z) 
M(a, 6 , z) 


< 


^z + 1 — 6 + 2 a + \/(z + 1 — 6 )^ + 4za^ 


Proof. Given H{z) 


M(a + 1 , 6 , z) 
M(a, 6 , z) 


, taking the derivative and using rules in 0 13.3] 


H'{z) 



6 + (z + 2 a + 1 — 6 )i/ (z) 


aH{zf) . 
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The characteristic roots of this Riccati equation are the roots ota + 1 — b+{z +2a+ 1 — 
b)X{z) — aA(z)^ = 0, which are 

(-z) = Y' (-2^ + 1 “ ^ + 2a ± ^/{z -\-l — bY + Aza^ . 

For z > 0 we have A+(z) > A_(z) and A^(z) > 0 if a > 0. These properties, together with 
the fact that H{z) = 1 + zjb + as z —^ 0 and H{z) ^ A+(z) ^ z/a as z —>■ +oo (see [9l 

13.7.1]) are enough to prove the theorem. 

Indeed, because H{z) is increasing for sufficiently small z (because H{z) = 1 + z/h + 0{z'^)) 
then A_(z) < H{z) < A+ for z small (because H'{z) < 0 if H{z) > A+(z) or H{z) > A+(z)); 
the same is true for large z because H{z) ^ z/a as z ^ +oo. But then necessarily A_(z) < 
H{z) < A+(z) for all z > 0 for the following reasons: 

1. Let Zc be the smallest value of z such that iL(zc) = A+(zc); we prove that such value can 
not exist. Indeed, if such Zc existed we would have H'{zc) = 0, H{zc) = A+(zc), H{z) < A+(z) 
for z < Zc and with A(^(z) > 0 for z > 0, which is not possible. Therefore H{z) < A+(z) for all 
z > 0 . 

2. For large enough z we have H{z) > A_(z) but then this holds for all z > 0. Observe 

that A_(z) is monotonic (increasing if a < 6 — 1, decreasing of a > 6 — 1, constant if a = 6 — 1). 
Then, if A_(z) is increasing or constant there can not exist Zc such that H{zc) = A_(zc) using 
the same arguments as in the item 1 of this proof. And if A_(z) is decreasing such number Zc 
can not exist either: if such a number existed then H{z~) < A_(z“) and < 0 but then 

this implies, because X-{x) is decreasing, that H'{z) < 0 for 0 < z < Zc, which is not true 
{H{z) is increasing close to z = 0 ). 

Then, because A_(z) < H{z) < A+(z) for all z > 0 we have H'{z) > 0 for all z > 0 □ 
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